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Abstract-Fluxes of cosmic ray particles with different LET values were measured on board the 

Cosmos-2044 biosatellite under various thicknesses of shielding by stacks of “ d 

plastic nuclear track detectors (mounted outside the satellite). The component compom of pamdes 

defected under shieldings of 0.1-2.5 g cm"’ is verified by comparing «P«nmen«l 

of model simulations of the fluxes of galactic cosmic ray particles and of radiation belt protons. 


1. INTRODUCTION 

The RADIATION fields in the Earth’s vicinity are gener- 
ated by the particles which are commonly grouped, 
according to some criteria (mainly their origin, com- 
position and energy) into certain components, 
namely, galactic and solar cosmic rays (GCR and 
SCR), radiation belt (RB) particles and anomalous 
fluxes (AF). To date, numerous experiments have 
been carried out on board Soviet and U.S. spacecraft, 
yielding copious data on the fluxes and spectra of 
protons and heavy ions in the Earth s vicinity (see, for 
example, Benton^and Parnell, 1988; Marenny er a/., , 
1987). Almost all of the data have been obtained 
using solid-state nuclear track detectors (SSNTD) 
(Marenny, 1987) which are notable for their essential 
characteristic of a detection threshold expressed by 
the smallest LET value at which detection is possible. 

In the present work, an attempt is made to analyze 
the components of SSNTD-detected particle fluxes 
under various thicknesses of shielding. The exper- 
imental data are supplemented with simulations cal- 
culated by models of the radiation environment. Only 
GCR and RB particles are included because other 
heavy cosmic ray components did not penetrate to 
greater shielding thicknesses than 0.05 g cm during 
the Cosmos-2044 flight. 


2. EXPERIMENT 

The Cosmos-2044 mission lasted for 13.8 days 
from 15 to 29 September 1989. The biosateilite orbit 


was elliptical, with a perigee and apogee of 216 and 
294 km, respectively, and an inclination of 82.3°. 
Four flat, lidded containers holding a variety of dosi- 
metric equipment were mounted outside the satellite. 

The experimental results presented here were ob- 
tained by processing three SSNTD stacks, one at 
IBMP (Moscow) and two at USF (San Francisco). 
The IBMP stack of 90 x 50 x 16 mm dimensions 
contained 20 Soviet-made KNC-type nitrocellulose 
detectors of 800 /nn thickness each. The detectors 
were etched in a 6N NaOH solution for 5 h at 50 C. 
The tracks were scanned with a stereo-microscope at 
80 x magnification. The penetrating (a cylinder or 
two cones) or single-cone (of at least 100 jim length) 
tracks were selected by counting, and correspond 
with the GCR particles of Z 5* 6 and LET Z 1600 
(MeVcm^g -1 (under the given etching conditions). 

The USF stacks were of 3 cm diameter and in- 
cluded CR-39 and Cronar polyester layers. The CR- 
39 SSNTDs were processed in a 6.25 N NaOH 
solution for 7 days at 50 °C. The bulk etch, was 
measured for each of the detectors. Pairs of detectors 
were reassembled in their flight orientations and the 
two adjacent inner surfaces were scanned with an 
optical microscope. This procedure permitted the 
particles to be separated into short-range (SR) 
(matching tracks appear on two inner surfaces only) 
and long-range (LR) (matching tracks appear on all 
four surfaces of the detector pairs) GCR particle 
tracks. The SR particles include the secondary SR 
particles from target nuclei within the plastic and the 
stopping primary GCR and trapped particles. All 
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protons were detected as SR particles because of their 
short registration ranges in CR-39. The GCR par- 
ticles include the primary GCRs and the LR Z 5=2 
secondaries which are mainly the GCR projectile 
fragments. The Z 2s 2 stopping GCR particles can 
also contribute to the SR tracks, but with a low 
probability. 

3. ANALYTICAL PROCEDURE 

3.1. Galactic cosmic rays 

The GCR particle fluxes outside the Earth’s mag- 
netosphere were calculated in terms of the model 
proposed by Suslov and Nymmik (1988) and de- 
scribed in detail by Suslov and Nymmik (1990). In 
the model, the energy spectra F,(E,t) of the 
1 Z < 28 GCR particles are inferred from the par- 
ticle rigidity spectra $,(£, t): 

F(E,t)dE = <P l (F,Oj^dR, (1) 


where the coefficients b, = 1.2 and d, =»> 0.034 describe 
the form of the spectrum at low energies 
(R < 0.55 GV) and were inferred from our data on 
the fluence of oxygen nuclei obtained in the Cosmos- 
2044 experiment (Marenny et al., 1990). 

3.2. Function of nuclei penetrating to low orbits 

The GCR particle energy spectra on the satellite 
orbit were found by calculating the penetrating func- 
tion on the basis of previous determinations of the 
boundaries for proton penetration into the magneto- 
sphere in a Cosmos-800 experiment (Biryukov et al., 

1984). . 

Figure 1 shows the plots of the penetration func- 
tion, calculated as indicated above and using the 
conventional techniques based on the International 
Reference Geomagnetic Field (IRGF) model Our 
approach has yielded higher fluences of particles 
penetrating to low orbits as compared with conven- 
tional techniques (Nymmik, 1991). 

3.3. Radiation belt protons 


where A i is the atomic number of a nucleus of species 
i; Rj is particle rigidity; t is time. 

The rigidity spectra are defined by the phenomeno- 
logical model 


W t) = 
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( 2 ) 


where D„ a, and •/, are the constants which character- 
ize the spectra of particles of certain species whose 
values were determined through the available set of 
experimental data on the GCR particle fluxes during 
the previous solar cycles. 

The data used to calculate the fluxes of all particles 
/-on be found in Standards for Galactic Cosmic Rays 
(1991); the constants used in calculating the energy 
spectra of the basic elements are presented in Table 1. 
The value of R, is defined by the solar activity level 
and by the delay of the CR particle flux variations 
relative to Wolf numbers. During the experiment, R* 
proved to be 0.88 GV. 

The power-law exponent in (2) is determined by the 
formula 
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Table 1. Values of constants used in calculating energy 


spectra 

Element 

Z, 

A, 

A 

a, 

7r 

H 

1 

1 

2.0E04 

3.0 

2.75 

He 

2 

4 

3.5E03 

3.0 

2.75 

c 

6 

12 

9.6E01 

3.1 

2.75 

o 

S 

16 

8.4E01 

3.0 

2.70 

Si 

14 

32 

1.2E01 

3.0 

2.65 

Fe 

26 

56 

9.2E00 

3.1 

2.60 


The RB proton fluxes on the Cosmos-2044 orbit 
were calculated using the AP-8 model (Sawyer and 
Vetter, 1976) for solar maximum. The calculations 
have shown that the orbit-integrated differential en- 
ergy spectrum of the 1-100 MeV protons can be 
described as 

F(E) = 8 x 10 3 £~ 7 , ( 4 ) 

where y = 2.03 at £ < 8 MeV and y - 2.72 at 
£ 5 = 8 MeV. 

3.4. Particle flux deep in matter 
The species i particle flux under shielding of thick- 
ness x is determined by the formula 

FJE'f) - F^LJLJi exp( - x /a,), (5) 
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where £'«£.- f* Ux) dx; L, and are the par- 
tide LET values above and under the shielding, 
respectively; A* is the species i parucle path fpr nuclear 
interaction. 

The particle energy spectra under the shielding 
were use d to find the particle LET spectra. 




where £,(£) is the energy dependence of species i 
particle LET value; j is a given flux component (GCR 
or RB). 

In conformity with the selection criteria for events 
in the USF stacks, the flux $,(>£) was broken into 
four groups (see Fig. 2): 


*,(>£)- t *S(>*)- (7) 

Jk-1 

If X, is the level of the upper surface of the upper 
detector in a pair of detectors, x 2 is the level of the 
lower surface of the lower detector in the pair, and x 
is the level between x, and x 2 , then the pictorial 
pattern for breaking into four groups of events is 
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where ® means that L > L„, O means that L < L„, 
and Y means that a particle is not present. The 



Fig. 2. The calculated flux vs detector depth. Curves 1-4 are 
particle fluxes calculated by formulas (8HH). respectively. 
S is the total panicle flux with L> L, at depth x. The curves 
numbered without a prime are for L >40 MeV cm g' ; 
those with a prime are for L > 1000 MeV cm g* . Curve 4 
is not present for the values of LET and plastic thickness 
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particle fluxes at L > L, at the level x may be 
described by 

*i(GCR) - Z<P,(L > L„ x, );(£(£) > EOh - *.)) 

* ( 8 ) 

4>i(SR) - £$„(£ 5* £., X, );(£(£) < £(x 2 - x,)) 

* (9) 

*i(SR) - !>*(£ < L » x,);(£(-R) < £(x : - x,)) 

* ( 10 ) 


<Pl (SR) - ( L<L » xi );(£(£) > ^(x 2 - *i»- 

9 (ID 

In the expressions (8 HU), the GCR group in- 
cludes the particles with LET values L > L, over a 
range great enough to penetrate both detector layers, 
from to x 7 . 


4. RESULTS AND DISCUSSION 

Figure 2 shows the (L > £»), plots $ * (GCR), 
<f>2(SR), #i(SR), and 4>1(SR) calculated for the 
experimental conditions under the shielding thick- 
nesses ranging from 0.075 to 2.5 cm H,0 at two 
threshold LET values (40 and 100 MeV cm' g - )• 
From the calculations it follows that, under the 
shielding exceeding 0.2 cm H 2 0, the fraction of the 
events in groups 2-4 (SR) is less than 1% of all the 
detected GCR particles, while the events of group 1 
include nearly the total fluence of detected CR par- 
ticles within the given component. Secondary par- 
ticles are not included in Fig. 2. 

Figure 3 shows the experimental data (plotted 
points) for (a) the particle fluences detected at seven 
depths at x > 0.2 cm H 2 0 in the USF stacks (GCR) 
at L > 100 and L > 500 MeV cm g"'; and (b) at two 
depths (0.24 and 1.53 cm) in the IMBP suck at 
l > 1600 MeV cm ! g~‘. The independence of the val- 
ues of GCR particle fluxes predicted by calculations 
for L>L, and at depths of 0.1-2.0gcm' : is 
confirmed by the dau. The measured particle flux at 
L > 100 MeVcnrg -1 is somewhat below that pre- 
dicted by the calculations, but probably within the 
overall accuracy expected from the comparison. 

Figure 3 also presents the dau on the events 
selected by the criteria of SR conditions (9H1 1) and 
shows the calculated dependence which follows from 
the RB proton spectrum. This dau (Curve 1) has 
been normalized to the SR particle flux at 
x = 0.25 cm H ; 0. To that end, the factor C in (4) had 
to be set equal to 1000, i.e. eight times less than the 
value given by the AP-8 model for solar maximum. 
The reduction in proton flux from the model predic- 
tion can be attributed to uncertainties deriving from 
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Fig. 3. Particle fluxes with L > L g vs detector depth. 
L > 100 meV cm g -1 : Curve 1 is the SR track density due to 
RB proton flux calculated by formula (4); Curve 2 is the 
track density due to GCR particie flux calculated by for- 
mulas (l)-{3); the black circles are experimental SR track 
densities; the circles with dots are experimental GCR track 
densities. L > 500 MeV cmg’ 1 : Curve 3 is the track density 
due to GCR panicle flux calculated by formulas (1)— (3); 
the star-like crosses are experimental GCR track densities; 
the skewed crosses are the SR track densities. 
L > 2 x 10 3 MeV cm g’ 1 : Curve 4 is the GCR track density 
calculated by formulas (1) — (3); the squares are the exper- 
imental GCR track densities. 

two major causes. There are substantial proton flux 
fluctuations observed at low altitudes near the lower 
RB border, where the density of residual atmosphere 
increases. The high solar activity level during the 
observations may have contributed to atmospheric 
expansion and a decrease in proton flux. Also, the 
pronounced directionality of the trapped proton 
fluxes is not taken into account by the AP-8 model or 
compensated for by the SSNTDs. The assumption of 
isotropic proton incidence which is made in correct- 
ing the experimental measurements for detector solid 
angle response could therefore have introduced some 
error into the measured proton fluxes. 

From Fig. 3 it follows that the SR particle flux at 
x >0.3 cm H 3 0 is considerably greater than would 
be expected from the primary RB protons. For the 
panicle flux with L > 100 MeV cm 2 g* 1 , the excess 
equals about 1 panicle (m 2 s sr)" 1 . Only a small pan 
of this can be attributed to H recoils produced by 
neutron collisions. The high energy (> 1 MeV) neu- 
tron flux was measured on Cosmos-2044 (Dudkin et 
a/., 1992) in the range of 275-500 r^j m" 2 s' 1 sr“ ! . The 
efficiency of CR-39 for high energy neutron flux 
detection has been reported as ~\0~ 4 proton tracks 


(neutron) -1 (Durrani and Bull, 1987). This would 
imply a neutron-induced component of -0.04 parti- 
cles m -2 s" 1 sr~ J , High energy protons also contrib- 
ute through elastic scattering of target H nuclei and 
this would seem to be the major source of the excess 
tracks seen at the larger shielding depths. 

The SR particles at L > 500 MeV cm 2 g" 1 cannot 
be protons, but are mainly secondary particles of 
higher charges produced by inelastic interactions of 
RB protons and neutrons with the detector nuclei 
(H, C, O). 
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